THEORY OF SITE-DISORDERED MAGNETS 



Th.M. NieuwenhuizenW and C.N. A. van DuW 2 ) 
™ University of Amsterdam, Van der Waals-Zeeman Instituut 
Valckenierstraat 65-67, 1018 XE Amsterdam, The Netherlands 
( 2 > Institute Lorentz for Theoretical Physics, Leiden University 
P.O.B. 9506, 2300 RA Leiden, The Netherlands 
(February 1, 2008) 

In realistic spinglasses, such as CuMn, AuFe and EuSrS, magnetic atoms are located at random 
positions. Their couplings are determined by their relative positions. For such systems a field theory 
is formulated. In certain limits it reduces to the Hopfield model, the Sherrington-Kirkpatrick model, 
and the Viana-Bray model. The model has a percolation transition, while for RKKY couplings the 
"concentration scaling" T g ~ c occurs. 

Within the Gaussian approximation the Ginzburg-Landau expansion is considered in the cluster- 
glass phase, that is to say, for not too small concentrations. Near special points, the prefactor of 
the cubic term, or the one of the replica-symmetry- breaking quartic term, may go through zero. 
Around such points new spin glass phases are found. 
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I. INTRODUCTION 

Prototypes of systems exhibiting a spin glass phase are the impure metals Cui_ c Mn c and Aui_ c Fe c , and the 
insulator Eui_ c Sr c S. [jjj Their common property is frustration of magnetic bonds, which arises from the combination 
of site disorder and the presence of both ferromagnetic and antiferromagnetic exchange constants. In metals the 
exchange is the long range, oscillating RKKY interaction; in an insulating SG the interaction is often antiferromagnetic 
for nearest neighbors and ferromagnetic for next nearest neighbors, while vanishing for more distant pairs. 

Edwards and Anderson proposed to describe the spinglass phase by spins with fully random couplings The 
mean field limit of this model for Ising spins is the Sherrington-Kirkpatrick (SK) model Q. The diluted version, 
where only a fraction of the random bonds is present, is called the Viana-Bray model jH). 

The mean field structure of the SG-phase was derived by Parisi and coworkers. There are infinitely many ther- 
modynamic states, and the order parameter is a function q(x) or P(q), which describes the probabilities of overlaps 
between pairs of states. Q 

The calculation of corrections to mean field theory has been very laborious. For a recent review, see Q|. 

The question what remains of mean field theory in low dimensions is still fiercely debated. Nevertheless, the best 
numerical simulations done so far are well analyzed within a mean field picture 

Mean field notions such as branching ratios and ultrametric organization of states have turned out useful in the 
interpretation of experiments. Indeed, noise experiments on mesoscopic metallic spin glasses have supported the mean 
field picture Q. 

We are thus at the stage where theoretical results are mainly gained from random bond systems, while experiments 
are mostly performed on random site systems. The understanding of the spin glass phase in random site systems has 
long been a challenge for the field. Some time ago one of the authors formulated a field theoretic approach for this 
purpose []To|| . In certain limit close analogies were established with known spin glass phases in well known models, 
such as the Sherrington-Kirkpatrick model, the Hopfield model, and the Viana-Bray model. 

The first purpose of the present work is to provide the details of that work. This will be done in section 2. In 
section 3 one of the results will be rederived via the Gaussian variational approach. In section 4 we consider a related 
model with combined bond- and site disorder. In section 5 we extend the approach of section 2 to vector spins. In 
section 6 we analyze the equations of section 2 for low T and find that replica symmetry must always be broken. 

More recently we have considered the Ginzburg-Landau theory of site-disordered spin glasses |ll| . It was argued 
that the prefactor of the cubic or most relevant quartic term may vanish at special points and then change sign. 
Around such points new types of spin glass phases occur. In section 7 we shall provide the details of that approach. 
In section 8 we close with a summary. 
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II. THERMODYNAMICS AS A SET OF COUPLED ORDER PARAMETERS 



We discuss a field theory for a random magnet with site-disorder. First we consider a system of Ising spins with 
Hamiltonian 

W(s) = — — J(r — r')s r s r 'C r c r i — H s r c r (1) 



r r' 



with sums over all N sites of a regular lattice in d dimensions, with periodic boundary conditions. The c r are random 
occupation numbers; c r — 1 indicates that at site r a spin is present, whereas c r — when this site is empty. We 
shall consider uncorrelated disorder with a fraction c and a total number cN of spins present, with < c < 1. 

r — r') is a translationally invariant pair coupling, a well known example being the RKKY coupling of 
metallic spin glasses, J(r) ~ cos(2kpr)/r 3 . However, we can consider more general situations. In a metallic SG there 
may be additional short-distance (anti-) ferromagnetic clustering For many insulating spin glasses J(r) takes the 
value Jnn > for nearest neighbors and J nnn < for next-nearest neighbors (nnn), while it is effectively zero for 
more distant pairs. 

In the following we shall focus on the situation where for large c ferromagnetic (FM) ordering occurs. For describing 
two-sublattice antiferromagnetic ordering one first has to make the replacement s r — > —s r on one of the sublattices. 
This will redefine the couplings J and introduce the staggered external field ±H. Since our spins are classical, this 
leads to a problem similar to the ferromagnetic one. 

In order to study the thermodynamics of this system at temperature T = 1/(3 we express the Boltzmann factor as 
a multiple integral 

exp{-f3H(s)} = J D<j)ex-p{-[3H{(j),s)} (2) 

where 

D(f> = det(/3J/27r + i0) 1/2 Y[ j d<t> r (3) 

r •* 

The integrations run from — \fioc to \fioc. To assure convergence of the Gaussian integrals, a regulator iO has been 
added to J; it will from now on not be written explicitely. The partition sum can thus be represented as trace over 
both discrete (s r ) and continuous (4> r ) degrees of freedom 

Z = J D<j>J2 e M-W(<t>,s)}, (4) 



with Hamiltonian 



rr' r 

where h = (3H. 

The quenched averaged free energy follows from log Z . The replica method is employed for studying this quantity. 
One thus calculates the quenched averages Z n for n = 1, 2, 3, • • •. The results have to be continued to the limit n — »• 0. 
This procedure introduces replicated field and spin variables </>", s" = s", for 1 < a < n. The average over the 
occupation variables c r can now be performed. Since spin sums are decoupled at each site, we may now write s™ as 
s„, which leads to 



with replicated Hamiltonian 



D0exp{-/3W„(0)}, (6) 



PHnW = \ it zZ ^J-M* ~ ciV$„. (7) 

a—l rr' 



where 
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= nlog2+ -^^]log ( 1 - c + ctr s exp^(<^ + h)s a J 

r \ a— 1 / 



(8) 



and where tr s stands for the normalized sum over s a — il; (p& — 1) * * * ? ^) 



tr * = 4 E = |r E ••• E ■ 

{s Q = ±l} Sl=±l s„=±l 

In order to proceed one has to make an assumption for the order parameter (s). In a mean field approach the 
magnetization per spin, M, will be proportional to the mean field value of m a — [(f> a ], where 

W^N-^^r (10) 

T" 

denotes the spatial average of an observable A r . The main step in the present work is to introduce, in the spirit of 
Edwards- Anderson, also space-independent composite order parameters such as [<j> a <i>p\c- This involves the second 
spatial cumulant 

Pa(3 = [^a^plc = [<Pa(t>p] ~ [4>a\ [^p] ■ (H) 

m a and p a p are the first terms in a cumulant expansion in powers of <p a . In Fourier space m a equals cj> a (0) and p a p — 
Ylk^o 4>ot{k)(j)p{— k). Ferromagnetic ordering implies a macroscopic occupation of the k = mode, </> Q (0) = 0(1), 
while a spin density wave has a macroscopic occupation of some k ^ 0. Spin glass ordering, however, occurs when 
Pa/3 > for a ^ [3 and is caused by a small occupation of all wavenumbers k =/= 0. 

In our approach the order parameters explicitly take into account the macroscopic contributions of combinations of 
Fourier modes, (e.g. M a = (f> a (k = 0)). The resulting Fourier sums (e.g. ^2 k ^ 4>a(k)4>/3(k)) "cannot bite", as their 
"teeth" have been removed. Likewise, at third order 

[4> a ^/34>-y]c = E ^a(ki)4>p(k 2 )<p 7 (~k 1 - k 2 ) (12) 

fc 1 ^0;fc 2 #0;fei#-fc 2 

only incorporates small contributions from all ki and k%, and no macroscopic (M a or q a j3) contribution. 

For a full description of the problem an infinity of spatial cumulants is needed. We shall, however, first see in how 
far low order approximations already lead to meaningful results. 



1. Second order cumulant expansion or Gaussian approximation 

As in the derivation of the Fokker-Planck equation, one may hope that the Gaussian approximation already leads 
to meaningful insights, and does not yet produce unphysical effects such as negative entropies. 

After expanding the logarithm of eq. (^) in powers of c/(l — c), we introduce in the Pth term the concentration 
factor 

^ = W^f ('=W,..) (13) 

and we have to replicate the spins I times. This leads to the "^-spin" 

i 

^ = E s « } - ( 14 ) 

The normalized trace over its values will be denoted by tr^ . When truncating the cumulant expansion beyond second 
order, the spatial average leads to 

exp^</> 



exp[exp 



- ex. u a\c 



exp 



(15) 



<>ri 
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The relation between p a p and the (fi's is imposed by inserting for each set a < (3 



d PaP H ^ dqape i N ^^-^\ (16) 

A Gaussian form in the <^>'s is then obtained, so they can be integrated out. This yields for the replicated free energy 
per spin (3F n = {-l/Nc)hxZ n : 

/^^E^+^/IIW 1 ^^ (17) 



a/3 



with 



OC 

$„ = n log 2 + ^ j e (tr< £) exp A (£) - l) . (18) 

Further, 

XW = y^(m a + /t)(7 a + ^ y^p a /3Cr Q ,(T / 3, (19) 

a a/3 

The fc-integral has arisen from (1/N) YlkjtO' wnere the exclusion of the k = term is due to the definition of the 
spatial cumulant. In the thermodynamic limit this exclusion becomes irrelevant. The ensemble averaged free energy 
per spin now follows as F — lim n _>o F n /n, with F n calculated at its saddlepoint. For c > \ the ^-sum in eq. ((lj) can 
defined by a Pade resummation. 

The magnetization per spin in the replicated system is 

T 

M a = -^—m a = < <r a > (20) 
cJ(0) 

OO 

= 5^trWo- a exD*W (21) 



i=i 



The physical interpretation of q follows from its mean field equation. This turns out to be very similar to the one 
introduced by Edwards- Anderson and Sherrington-Kirkpatrick 

qafi = (v a cr(,) = 2^7£tr^V Qt j^expA w = ^ le ^ — - . (22) 



i=i 



tr^expA^ 



The last equality holds since the denominator equals unity in the replica limit n — > 0. For c — > only the £ = 1 term 
(s a s/3) survives and one exactly recovers the two state overlap introduced by Edwards- Anderson and Sherrington- 
Kirkpatrick. Variation with respect to q a p yields 

„-f"(-jM-) . ,23, 



(2ir) d I 1 - c/3J(k)q 



a/3 



This form exhibits clustering effects. Indeed, decomposing q a p — qdS a f3 + <jaf3 with qd — q aa and q aa = 0, we may 
also write for a ^ (3 

f R t (U\ ( 1 ^ 



\i- c pj oS {k) q> 



a/3 



c^J^(fc) -4— - (« ± P), (24) 



(27r) d " cnw I 1 - c/3J eS (k)q 



a/3 
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with the effective coupling 



J eff (fc) = J -M . (25) 

Expansion in powers of c indeed shows that it is a multi-spin effect. Maxima of J are more pronounced in J e g. This 
describes the formation of finite clusters that are precursors of the would-be spin density wave or (anti-) ferromagnetic 
phase transition. In metallic spin glasses incomplete spin density waves have indeed been observed by neutron 
scattering. jl2| It is well known that both in metallic and in insulating spin glasses there appear ferromagnetic clusters 
near the ferromagnetic transition line; these clusters (sometimes called "fat spins") act as a quite rigid effective spins, 
that may contain up to 2000 magnetic atoms, making quite visible mov es. They are responsible for the dynamics in 
the cluster glass phase jjj . We shall return to this point in Section VII . 



(At low temperature this expression for J c ff becomes singular, though eq. ( |23| ) is well behaved. It is then more 

appropriate to interpret J/(l — c/3J(qd — Qea)) as the effective coupling, where qEA = max q a ^p is the Edwards- 
Anderson order parameter). 

2. Replica symmetry 

To obtain an idea of the content of previous expressions, we consider them in the replica-symmetric sector. The 
free energy then reads 

flP T 2 1 r d d fc /, „ Qinw „ c(3J(k) 
2cJ(0) 2c 
1 

+ 2 



pq) - ^ 7l J g{x)dx\n J g(y) dy2 e cosr/(/i + m + x^/p + y^Jp d - p), (26) 



where 



e~ x I 2 



is the Gaussian weight. The saddle point equations read 

d d k cl3 2 J 2 {k) 



P 



Pd=P + 



(2TT) d (l- C f3J(k)(q d -q)) 2 ' 
d d k (3J{k) 
(^) d l-c(3J{k){q d ~ q y 

J g(y) dy cosh^ -0 tanh ip 



Q = ^2li / g{x)dx£ 



J g(y) cosh v 



V- f , /g^dj/cosh^^ + ^-^tanh 2 ^;) 

Qd = 2,11 I g(x) dx -j- , (28) 

t J J g[y)dycosh ip 



where 



ip — h + m + x\/n + y\Jpd — P- (29) 
These expressions will be used later on for small T, where the cosh's essentially become exponentials. 

3. Relation to the Hopfield model and the Sherrington- Kirkpatrick model 

Despite of the simplifications made, the above expressions are very rich. In some limits of long range couplings 
they become exact. 
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Consider the situation where 



J(k) = Jo k < \k\ < ki (30) 

and zero outside this shell. This corresponds to a long range spatial coupling 

f d d k sin(fcr) 
J{r) = Jo ttt-t-t — : • (31) 



fcn <fc< fcl P*) d kr 

We denote the volume in phase space of this shell by A = d d k/ (27r) rf . Consider the combined limit c — > 0, k± — » ko, 
such that a = A/c remains fixed. Inserting this in eq. p7[) we obtain for its spin glass content (m = 0) 



a otf) 

In order to understand its meaning, let us consider the case of replica symmetry. We get in the limit n — ► 

0F - | {.o g (l - « - ,)) - 1 _ < g^_, ) } + - .) 

'fe =e -x 2 /2 log2cosh2 ,^ (33) 



with p and g obeying 



/2tt 



P 



a(c/3J ) 2 g 



(1 - c/3J (l - <?)) 2 
da; 



<1 



'2tt 



e~ x /2 tanh 2 Xy/p. (34) 



After a rescaling (T — > cJqT, p — > ar/T 2 ) these are exactly the replica-symmetric free energy and its saddle point 
equations for the spin glass phase in the Hopfield model of a neural network with a N patterns jl3| . More generally, 
eq. ( [32] ) is the free energy including replica symmetry breaking. In the limit A — > 0, c — > 0, with fixed A/c = a one 
considers long range interactions between widely separated spins, that effectively lead to random interactions. It has 
been checked that loop-corrections do not contribute in the limit. So (^) is an exact mean field theory in this long 
range limit. Note that this may occur in any fixed dimension (d = 1, 2, 3, • • •). 

The SK-model arises in the subsequent limit a — > 00. To show this we must rescale T — > TcJo^fa. Then the first 
relation in eq. ( |34| ) becomes linear, and eq. (|32|) and ( |33"| ) become linear+quadratic in q a /3, so that p a p and q a p are 
linearly related. After elimination the p's, we recover expressions that in the limit n — ► are equivalent to the SK 
model. Hereto we also use that for small n our expression — 1 + tr s exp(p Q ^s Q S/3/2) is proportional to n and to this 
order equivalent to the SK result logtr s exp(p ct ^s Q S/3/2). 

The present random site problem was also studied, for the case of Heisenberg spins, by Bray and Moore |Q. These 
authors do not introduce conjugated spin variables but make a cumulant expansion of the Boltzmann factor up to 
second order in the occupation factors Cj. Since the latter are zero or one, there is, however, no reason to stop at 
second order. Therefore the results differ from ours and are incorrect. For instance, according to ref. fl4|| the spin 
glass transition temperature scales as T g ~ c for small c both for short range and RKKY couplings. We find, however, 



from eq. (34) that T g — Jr ~ *Jc. Below we shall derive T g ~ c for RKKY couplings via a more intricate analysis, 



that takes into account more contributions. 



4- Full cumulant expansion: the direct way 

The above results do not exhibit a percolation transition. Indeed, T sg ~ y/£ does not vanish identically for small 
enough concentrations. This is related to the fact that the percolation transition of systems with z — 00 occurs at 
c = 0. Also for RKKY couplings the present prediction overestimates T g . Here one expects the well established 
"concentration scaling" T sg ~ c. jjj 

We therefore consider the full cumulant expansion of eq. (JsJ) . We introduce order parameters 

Pax-a h = [4><xx ■••<Aa fc ]c (35) 



G 



as well as conjugated parameters q ai - a k - They are symmetric in their replica indices. The replicated free energy may 
be written as 

(3F n = +pq- $ n (p), (36) 

where 



\ arr' fc=l ar-Qfc 



(37) 



and 

w = X]w ■••«> '/■• • ( 38 ) 

fc=l ' a r -a t 

<&„ is given by eq. (fq) with 



2| 3 , 

a a/3 a/37 

The free energy is now determined by the infinite set of conpled equations 



d 

uq ai ... ak 
d 

q ai -a k = k\- * n (p). (40) 

OPat—ati, 



5. A resummation 

The above result is not very useful at low temperatures, where the p's become very large. However, it is possible 
to resum the contributions. This can be performed in the following elegant way. In the £'th term of <!>„ there are spin 
variables s«\ (J — 1, ■ • • ,£), which we now label as s a with a = (a,j). Though they do not depend on j, we also 
denote <f> a as <p a . We start from the definition 



[expj^ 4> a s a ] = exp [exp^ 4> a s a ]c- - 1 

a \ a / 



(41) 



The obvious relation 

[A + A a <j> a r + A ab tf<t> b r + A abc rA h A L r H }c = A + Aain + A ab [cj> a ^] c + A abc [<p a "cf> 1 '<p% + ■■■ (42) 

(summation convention employed) defines cumulants of any expandable local function of It is a generalization of 
the definition 

[e*° Sa ] c = 1 + S a [(j) a ] c + -S a S b [(t) a <j) b ] c + -S a S b S c [(j)a4> b (t>c]c + 7-rS a S b S c S d [(j) a (j) b (j) c (j)d]c H • (43) 

2 o 24 

We employ the identity 

exp(J^ (f> a s a ) = TJcosh0 a (l + s a tanh0 a ) = JJc„(l + s a t a ). (44) 

a a a 

Using the linear property of series of cumulants, eq. (|42|), this product brings 

[exp^ 4> a S a }c = [J]^C e ] c + ^ S Q [t a rjc e ] c + - ^ S a S b [t a t b ]~[ c e ] c 
a e a e a^b e 

~F X! s aS b S c [t a t b t c Y\_C e ]c+ -^J S a S b S c S d [tat b t c tdY\_ c e]c H ■ (45) 



a^b^c^a e a,b,c,d; different 
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That the subtracted terms indeed work out this way can be checked for low orders in the small <p a expansion. Up to 
quartic order the right hand side gives 

[i +\€ + \<t>l<i>l - ^&\c + M«(i - \<& + \<ft)]c 



+ ■^[SaSbS c Sd4>a4>b4>c4>d ~ 6s a Sb(f)a4>b<t>l + i^a^b + 8s Q Sfc0a0fc ~ 60*]e, (46) 

where now the sums over repeated indices are unrestricted. This expression indeed coincides with eq. (|4 
For every integer 1 < k < n and t > 1 we introduce symmetric order parameters 

Pan— = [ c l ' ' ' c n^OLi ' ' ' t ak ] c , (47) 



with the restriction that the total number of replica indices that occur more than once, does not exceed £. (For 

. Further 

pW = [el...ci] c -l. (48) 



example, p^l^ is needed, but is not). Further 



Introducing the spin variable 



2| / jfaO"^"" ' gj 

a^b a^b^c^a 



we find that the i'th term of involves an expression tts exp{p^ + X^}. Next we introduce the conjugate order 
parameters q£\ q^, • • •, through relations of the type (|l6|), 

1 = J dp ( al- ak J 'j^i dq °^- a » eX P^^ C 9ai---a fe ([*ai ' ' " II ~ ( 50 ) 

Let us first consider the saddle point equations for p^> and q^>. It is readily seen that p^' — 0(n) and that 
= 1 + 0(n). This implies that p^ can be omitted and that the q^ terms can be summed. These manipulations 
lead us to a replicated free energy of the form , 

0F n = -Hf n (q) + (p,q)-$ n (p), (51) 

with 3> n given by eqs. (|l8| ) and (^9|). Further, 



oo oo 



fc! 

-1 fc— 1 ar-'Qfc 



and 



* n = -^log y^exp ( ~ <P? J rM?> + N 



log I 1 - c + ctr s exp ^ <f> 



cN ^ 7£ ^ — ^ gai-a, [Cj • • • C e n t ai ■■■t ak 
£=1 k=l ' ct!---a k 



(53) 



In this expression the logarithm sums the q^ — > 1 contributions. The difficult part of the problem is still the 
integration over the 4> fields. One might even think that nothing has been gained in the present representation. 
However, the dangerous order parameters have been subtracted explicitely. In the remainder of the paper we shall 
mainly consider situations where has a relatively simple form. Nevertheless, we shall be able to draw interesting 
conclusions. 

In general, the (f> integrals cannot be performed exactly. If the exponent is expanded up to second order in (j>, one 
recovers eq. ([17]). 
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6. Small concentrations and the relation with the Viana-Bray model 



We now consider the limit of small c at fixed T. Here only £ = 1 contributes, so that all q's and p's with £ > 1 can 
be omitted. The remaining ones have all replica indices different from each other. ^> n can be approximated to order 
q 2 . The integrals can be performed. To see how this goes, consider eq. (^3|). First, since c is small, we may omit the 
term iV[log- • -] c in the exponent. When expanding in powers of q we have bilinear terms of the type 

c 2 q ai ---a k q a ' 1 --- a ' k , /-D0exp(-^ ^ <j% J~ r )(j)" )N[c 1 ■ ■ ■ c n t ai ■ ■ •iajc^ci • ■■c n t ol ' x ■ ■ -t a i^] c . (54) 

arr' 

Neglecting, for the moment, the difference between cumulants and ordinary moments, the integral can be written as 

f\ )T"2 s,s " arr' 7 

= ^2 ^2 s ai ■■ ■ s ak s a/i ■ ■ ■ s a > k/ exp(/3J(ri - r 2 ) ^ s 7 s 7 ). (55) 

ri,r2 s,s 7 

For each fixed 7, the sum over s 7 and s 7 can have preexponential factor 1, s 7 , s 7 , and s 7 s 7 . When performing the 
sums, it is readily see that only the first or the last type of terms survive; this explains also why terms linear in q 
vanish identically. This implies that only diagonal terms in the q ai ---a k remain, with prefactors that are spatial sums 
of cosh 71 /3J(ri — r 2 ) tanh fe /3J(n — r 2 ). Investigating the subtraction terms that occur in the spatial cumulants we 
observe that replacing the cumulants by ordinary moments leads at this point only to corrections of order 1/jV, that 
are negligible in the thermodynamic limit. The same effect led to the exclusion of the k = term in (Jf7|). This 
supports previous finding that the role of the order parameters is only to eliminate the macroscopic parts of spatial 
averages; the remainder being the well-behaved cumulants. For small n we find for the replicated free energy 

&Fn = § \ -riTQ + n^Ml + Ti Y,q 2 af3 +T 3 + \ -nlog2 + l-tr s expX, (56) 



2 

with 

T() 

and 



ct</3 a</3<7 



^2 log cosh /3J(r) r fc = tanh fc (3 J{r) (k > 1) (57) 



X = y^^CTlMg +h)s a + CT 2 ^2 1af3S a Sp + CT 3 ^ (58) 
a a<fS a</3<7 

When couplings are (mainly) ferromagnetic, a ferromagnetic transition sets in at temperature such that cri = 1. If 
the J(r) are partly positive and negative, there are cancellations in this sum. Then a spinglass phase may occur at a 
higher temperature. The spinglass phase (q a p > 0) sets in at a temperature T g where ct 2 = 1. If only couplings to 
z neighbors are different from zero, a SG phase can only occur beyond the percolation threshold c p — l/z. c/c p > 1 
expresses that, on the average, each spin should interact with more than one other spin. For a fee lattice in d = 3 
this yields c p = 1/12 w 8.3%, to be compared with the simulated value c p w 13.6% pl|. For the RKKY-coupling in 
three dimensions, the condition ct 2 = 1 indeed yields the "concentration scaling" Tq ~ c for small c. 

Some years ago the SK model with diluted random bonds was introduced by Viana and Bray ^|; it is commonly 
called the Viana-Bray model. These authors also analyzed various aspects of the phase diagram. In an unpublished 
work one of the present authors fl6|| derived the same model, without making any analysis of the phase diagram. In 
the diluted model one assumes that the couplings J rr * in eq. ([j]) vanish with probability 1 — p/N, and are drawn 
independently from a normalized distribution r( J) with probability p/N. In the thermodynamic limit TV — * 00, the 
parameter p is kept fixed. There is a close analogy with the Flory model of coagulation Q The parameter p can 
be interpreted as a dimensionless time variable, at which the coagulation has been interrupted. For p < 1 only finite 
clusters exist, while for p > 1 an infinite cluster has appeared. On this infinite cluster a spontaneous phase transition 
can occur, to a (anti-) ferromagnetic or spinglass phase. 

It would be interesting to consider p really as a slow time variable, and to consider the combination of very slow 
cluster growth and spin glass dynamics. In this dynamical model the random couplings would be fixed, but as time 
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progresses more and more of them would become active. This idea is inspired on a model, introduced by Coolen, 
Penney, and Sherrington, where all couplings are active, but their strength changes slowly in time |l9f| . 

The quenched free energy of the Viana-Bray model is derived along the following lines. The averaged replicated 
partition sum equals 



(59) 



Z '~E II fl-| + |/r(J)«l/Bp^J^)} 

{s°}<rr'> { a ) 

~E ex p|^E(- 1 + / K^)dJcosh"(/3J)J](l + «,tanh/3J)^|. 

Writing out the product over a and decoupling the terms quadratic in products of spins at the same site, one defines 
to = J r(J)c£7 log cosh /3J; t k = J r{J)dJ tanh fc (3 J (k = 1, 2, 3, • • •)• (60) 

In the limit n — > the replicated free energy takes the form 

l3Fn= 2 -^o + ^E^ + ^E^ 1 ' 3 E + -nlog2-logtr s expX, (61) 

[ a a<(3 Q</3<7 J 

with X given by eq. (p8[L after the replacement CTk — > pifc. For a discussion of the phase diagram of this model, see 
e.g. refs. §, @, 

For short range couplings in the random site problem the two results are thus closely related. Indeed, since 
again — 1 + tr s expX is equivalent to logtr s expX, the only real difference lies in the meaning of the parameters 
cTfc = (c/c p )(rfe/ z) = ptk- In the diluted random bond model Tk/ z — tk is defined as the average of tanh fc (3 J over the 
distribution of the independent couplings J. In the random site model with sure couplings, this average is taken over 
space. 

III. GAUSSIAN VARIATIONAL APPROACH 

Since disordered systems are so hard to solve, it is becoming more and more popular to investigate variational 
approaches p3| , [p4| , pq] , p6| . In the present problem one starts from eqs. (Q) and (JsJ) , and chooses the variational 
Hamiltonian 

^ = (0? - m a ){G- l }f{$ - mp) (62) 

The variational propagator (4>f<pj) vav = Gfj* = G a p{ri — Tj) is translationally invariant and has Fourier transform 
G a p{k). The variational free energy per spin becomes 

0F = — (- lnZ var + (H n - H var ) var ) 
Nc 

V- rn 2 a 1 ^ f d d k r , G(fc) , 

+ J2 7/(1 " trWe^= (,l+ma)CTQ + ^ ^ ^^=0)^ (63) 



=1 



The variational equation 8F / 8G a p{k) = can be expressed in terms of q a p using (|22|). This yields 

(3J(k) 
1 - c/3J(% 



a/3 



Its fc-integral can be identified with eq. (E3| 
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^ = G ^ = ») = /^( r ^) o; (65) 

Due to this, the £-sum in eq. (^) can be identified with the expressions (|l8|), (|l9|). Therefore the variational equations 
of the second order cumulant expansion |L0| and the Gaussian variational approximation coincide. It is then a small 
exercise to check that the saddle point value of the free energy also coincides. This does not mean that both approaches 
are the same. The Gaussian variational Ansatz has looked in a larger space, with variational parameters G a p(r), 
but found the same physics as the second order cumulant expansion with its space-independent parameters p a p and 
q a p. The reason hereto is simply that both approximations are Gaussian. The equivalence of the two approaches was 
realized by one of us in fall 1995. 

One can go to the continuum limit. Hereto one reinserts the lattice constant a, and takes the combined limit c — > 0, 
a — > 0, such that p = c/a d is fixed. The main effect is an overall factor a d (since in the limit F/a d becomes the finite 
free energy per unit volume) and that only the t = 1 term survives in the sum (because c->0). This continuum result 
was discovered independently by Dean and Lancaster |27j], who started in the continuum immediately. (We should 
point out that a proper definition of their path integral would reintroduce a lattice.) They introduce a grand-canonical 
description of disorder, which is physically the same as the ordinary disorder average, but for their purpose a bit more 
convenient. 



We have criticism to this work. The claim of 27 to present a new field theoretic approach is unjust. Their field 
theory in replica space (paper I, eq. (4), paper II eq. (2.5)) is just the continuum limit of the above eqs. (g), (0), (||), 



represented long ago as eqs. (2) and (3) of ref . [|10|. The authors of 
the different approach, the resulting saddle point equations and sadd 
coincide with the ones of the second order cumulant expansion of |10 



27[ fully refrained from explaining that, despite 
c point free energy are not new either |28[ , but 
I, and also discussed above. 



IV. A SITE-DISORDERED SPIN GLASS MODEL WITH PAIR OVERLAPS ONLY 

In realistic spin glasses the signs of the interactions oscillate due to the strong distance-dependence for pairs of 
spins present. We have seen that this leads to global multi-spin order parameters, a somewhat unfamiliar theoretical 
framework. On the other hand, it is custommary to consider models where the signs are uncorrelated random variables. 
We shall extend that approach to the site-diluted case, and find that only local pair overlaps occur. 

Consider the couplings 

J rr , — J(r — r')£ r £ r > (66) 

with independent quenched random numbers £ r . Numerically the simplest choice is £ r = ±1 with equal probabilities. 
We shall take them, however, Gaussian with zero average and unit variance. The function J(r) can now be taken non- 
negative. For J(r) ~ 1/r 3 this form replaces the deterministic signs of the RKKY interaction by random signs. We 
proceed along the lines of section 2, and introduce fields <fi" . The disorder average over the random spin configurations 
and the random £'s will now yield 



exp ( ^ S a C r (4>r& +h)\ = 1 - C + C CXp j ^ ^ <t>r <Pr s a$ + h^S a \ 
\ a / \ a/3 a J 



(67) 



We can introduce the composite local field p a p(r) = </>"<^f as new variable by inserting 5 functions at each lattice site. 
In their plane wave representations there appear only Gaussian <f> integrals. Integrating them out one is left with a 
field theory for the space-dependent fields p a p(r), q a p(r). Its Hamiltonian reads 

pHn = 0H(m(r),p(r),q(r)) (68) 

r 

with 



In [ 1 - c + ctr s exp ( - ^2p a ps a sp + /i^]s a j 



(69) 
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Note that the log-term is non-local. In the mean field approximation, the functions p and q are space- independent. 
The resulting expression is then very close to the one in the second order cumulant expansion of section 2. Now the 
spin sums are simpler, since the logarithm n eeds not be expanded. Consequently, the sign changes in the coefficients 



w and j/i, 2/2, 2/3, to be discussed in section VII, do not occur in the present model 



As in section 2, eq. rt29), the effective coupling is 



Jeff(fc) - ^ , (70) 

l-c(3J(k)(l-q EA y 



which again exhibits clustering effects. 



V. VECTOR SPINS 



Experimentally, spinglasses usually consist of Heisenbcrg spins with weak or strong anisotropy. It is therefore useful 
to extend our formalism to m-component spins. We consider the Hamiltonian 



H(a) = ~ E ■CW'V'V - E H r^c r (71) 

fivrr' fir 

where |S| = 1. The translationally invariant coupling J may contain effects of anisotropy, such as dipolar anisotropy. 



7. Second order cumulant expansion 
Proceeding as above we can carry out the calculations in the second order cumulant expansion. We find 

^ = Y c £ + Yc I E^ 1 " «%)W W + \ E pX ~ ( ?2 ) 

The expression for takes the form 

oo 

$„ = n log fi m + le{ - 1 + tr^ exp } . (73) 

Here fi m is the area of a unit sphere in m dimensions (fii = 2, 0,2 = 2n; Q3 = Ait). DS = r2~™c?Si • • -dS n , is the 
angular integration measure for one replicated spin and tr^ = Y\j=i I DS^ is the integration over all replicated 
spins. Finally, 

XW = E« + £ Sjp" + \ J2 jft jz S<p"Sjf» (74) 

ap, j=l afSpv jj' = l 

In the SK-type limit discussed above, one again obtains and expression quadratic in the c/'s: 

^ = | E + 15 E E E + \ E pX - ( 75 ) 

a.fivr r a/3 pvper a(3pu 

For isotropic couplings in the presence of an external field, Gabay and Toulouse ^| already showed the occurrence 
of an irreversibility line where transverse components of the spins freeze. 

8. The resummed cumulant expansion 



Here our aim is to derive the equivalent for vector spins of the resummed cumulant expansion of section [I 5 . This 
will yield, in particular, a prediction for the transition temperature at low concentrations. 

The first step in the derivation is the vector analog of the cumulant generating function eq. (Ell): 
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[cxp £ = exp ( [exp £ #S£] C - 1 ) 

ail \ a/i / 



(76) 



where = 1, 2, • • • , m and where a = (a, j) with a = 1, 2, • • • , n and j = 1, 2, • • • , I. The natural generalization of eq. 
( p4| ) is to introduce 



exp^T 0£5£ = C a + C a T a (S a ) = P a exp^ <(£S£ + Q a exp^ 



(77) 



where 



C a = -^~ f dS a exp £ = ^ exp £ 

with P a = fJ^ 1 J <iS a being a projector, Q a = I — P a being its complement, and 

rp (a \ — cx P^m<^ , 

J-a(^a) — -~ J- 



(78) 



In analogy with section II 5 we define order parameter functions 

Pai—ak (S Q1 , • • • , S afc ) = \C[ ■ ■ ■ CnT ai (S ai ) • • • T ak (S ak/ 



(79) 



(80) 



where the total number of repeated replica indices should not exceed I. Note that this qua ntity depends only on k 
spin variables, generalizing a similar property of the combination p ai ... ak s ai ■ ■ ■ s ak of section II 5. The expression for 
$„ takes the form ( 173|) with 



xW=^{pW(S a ) + h.S a } + i^S(S a) S 6 ) + i £ pW y (S.,S 6) S e ) + ... 



3! 



Going through similar steps as before we derive the replicated free energy 

f3F n = -*„(<?) +pq-®n(p), 

where 



(81) 



(82) 



* n =-Llog /^exp(-|^<(r){J-Var(r)+7V log ll - c + c j DSop^^ 

(oo oo 1 \ 

c7V E 7< E M E P l • • • ^ (Si - • • • > • • • Ci n T "i ( S l) • • • ( S ^)] c • 

£=1 fc=l ' ai-o t / 



and where 



oo x ! 

M = E 7£ Ew E p i'-- p fePai---a fc (Si,---,S fe )gi £ 1 ) ... afc (Si,---,S fc ). 



fc! 

£=1 fe=l oi'"a t 



(83) 



(84) 



In the limit of strong anisotropy, where the vector spins are Ising-like and point along the ±z axis, these equations 
reduce to the results of section [I 5. Indeed, it can be seen that qa}... ak (Si, • • • , S&) — > {qfl— ah } isingSf ' ' ' m this 



limit. 



9. TTie spin g/ass transition at low concentrations 

The spin glass transition occurs for small p a g and q a L We now restrict ourselves to the case of small c, where 
only the 1 = 1 terms contribute. The transition temperature is found by considering a replica symmetric solution 
<Zq/3(Si, S2) = <Zsg(Si, S2). By expanding the mean field equations oip a p and q a p it is seen that the condition for the 
critical point has the form of an integral equation 
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P 3 P 4 9sg(S 3 , S 4 )gi0 2 Q3Q4 exp feuAstPj^sg + sspjgsz}) 

q SG (S 1 ,S 2 )=cJ2 V — 2 L , (85) 

{P 5 P6exp(E^^/3J«)} 

where J(f" — \{Jq" + J"o)- Note that this equation has the form eigenvalue equation, Kernel * qsa — QSG- This 
prediction is useful for long range spin glasses with a small concentration of spins, as typically occurs in experiments 
on metallic spinglasses and also in the Monte Carlo experiment of ref. ||(J . 

The fact that an integral equation occurs is not unexpected, since this is already known to happen for a linear 
chain. Because of this reason, the analysis of the SG phase for very diluted vector spins is more cumbersome than for 
Ising spins. 

For isotropic couplings, jv(r) = S^J(r) the relevant solution of eq. ( p5| ) has the form ?sg(Si,S2) = Si • S 2 . 
Indeed, the QiS may be replaced by unity because of symmetry and the remaining integrals over S3 and S4 are 
proportional to Si and S2, respectively. This leads to the condition 

JdSS z ex P ((3J(r)S z )\ 2 _ ^ J I f (/3J(r)) 



44 /dSexp09J(r)S«) J ^\l^i{PJ{r))\ ' [ ' 

where I„ is the modified Bessel function of inde x v. It is the generalization of the relation ct 2 = c^ r tanh 2 f3J(r) = 1 



for Ising spins (m =1), derived in section [II q . For Heisenberg spins (to = 3) eq. (86) becomes 

f 1 ^ 2 

c^|cotanh/3J(r)-^^ 



= !• (87) 



As it involves the classical Langevin function for our classical spins, one expects it to contain to the quantum Langevin 
function for quantum spins. Eq. ( |S7| ) is quite natural, as the same mean field form, but without taking the square, is 
known to occur for the ferromagnetic transition. 

Eq. ( p7| ) again leads to concentration scaling Tq ~ c for RKKY couplings. We have tested it for the situation 
considered by Matsubara and Iguchi. ^TJ| These authors consider a fee lattice with lattice constant ao and with c = 5% 
spins present. The RKKY coupling J(r) = Jq cos(2kFr)(ao/r) 3 involves the parameter kp = 4.91/ao; J(r) is set to 
zero for r > 3a . The simulations yield a spin glass transition at T g /J = 0.068 ± 0.008. We find from the above 
relation that T g / Jq = 0.306. It shows that c = 5% is not a small concentration, and that higher order corrections in 
c are important. 



VI. SPIN GLASS BEHAVIOR AT LOW TEMPERATURES 

In order to see whether spin glass behavior in the form of breaking of replica symmetry is likely to occur, we 
consider the simplest generalization of the SK- model. It is given by eq. (|l7|), (p!s|), and |l9|). The model represents 
the site-diluted spin system in a Gaussian approximation. 

Consider the low temperature behavior of the replica symmetric solution M a = M, q aa — q e i, q a ^i3 — <7- For small 

T we observe the "groundstate dominance" ss) — si, 1 ' for all j = !,-■■,£. This implies that repeated spin sums 
reduce to single spin sums. It leads to a free energy 

T , 1 f d d k [ , ~ . . cBJ(k)q } 

PF = s to 2 + — / — — r <^ In 1 - c(3J(k)(q d - q)) " V n \ 

2cJ(0) 2c J (2nY\ 1 Kmd 1 - c0J(k)(q d - q) j 

+ niPdQd ~ pq) - ^ ll J g(x)dx\n2cosh£(h + m + Xy/p), (88) 



where g(x) = exp(— x 2 /2)/\/2tt is the Gaussian weight. It can be checked that q d = 1 — c. The other mean field 
equations are 



M = J 11 = [ q(x)dx'S~ r -ff£t&nh£(h + to + x x fp], 
cJ(Q) J V 

q = lg(x)dx'£ 'It? t«»h 2 l(h + m + xjp), 

J n 
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d d k c/3 2 J 2 (k) 

(2T) d (l-c/3J(A;)(gd-9)) 2 ' 
f d d k /3J(k) 

Pd=P+ 77m- 7TT77T7 ( 89 ) 



(2ir) d l-cf3J(k)(q d -q) 
For small T it holds that p ~ T~ 2 is large. One therefore finds 

q = Qd - \ - z2^jn= = id - J — . (90) 

Note that this result also holds for c > 1/2, where the £-sum is formally divergent. 

According to the analysis by de Almeida and Thouless f3l| , the instability of the replica-symmetric solution is 
caused by a negative eigenvalue of the fluctuation matrix. This mode has been termed "replicon" . |5^] We have 
proposed the name "ergodon" for it, since this name relates to the underlying physical mechanism of ergodicity 
breaking, rather than to one of the mathematical ways to study it. |33) The "mass" of the "ergodon" or "replicon" 
can be derived along the lines of AT. One expands the free energy up to second order in 6p a p, Sq a p, imposing the 
condition ^2 a 5p a p — Yl a = ^ This finally leads to a 2 times 2 matrix, which has one large eigenvalue (not of 
interest to us) and an eigenvalue that may become zero and negative. The latter eigenvalue has in case of ground 
state dominance the approximate form 



Using q — > q d = 1 — c, this leads to 



A(T) w 1 — (1 — 2c)\ —. (92) 



As p ~ T -2 , it follows that A ~ —(1 — 2c) /T, showing the expected AT instability for < c < 1/2. However, for 
c > i the replica symmetric mode is seen to stabilize at low T . 

To see whether the replica symmetric state has correct physics for c > 1/2, we consider the zero-point entropy. We 
take m = h = 0. To order T° the expression (B8J) now becomes 



p F = I [ 4% lln(l - c(3J(k)( qd -q)) \ - (93) 

2c J (27r)«n 1 1 ' W q " l-cpJ{k){q d -q) ( 1 




It follows that the entropy has the form 

1 f I 1/1 flfaw ^ c0j(k)(q d -q) \ ^ n 

$ = 7T 77T\d S - ln ( : _ c P J K k )Kld - q)) - -st—- ) < 0. (94) 

2c ^ ( 27r ) [ l-c/3J(k)(q d -q) J 

For T — > this goes to a finite but negative value. The same behavior occurs in the replica symmetric solution of the 
SK model, which follows as a limiting case of the present model. 

Thus replica symmetry must be broken for all c at low T . There is a glassy phase, possibly coexisting with (anti-) 
ferromagnetism. This is somewhat surprising when the replica symmetric solution is stable, but similar behavior was 
found in a neural network. |34| The present calculation suggests that infinite order replica symmetry breaking occurs 
at low T. Due to the "ground state dominance" , it could be analyzed along the lines familiar from the SK model. 

VII. GINZBURG-LANDAU EXPANSION NEAR THE SPINGLASS TRANSITION 

We use a Ginzburg-Landau expansion of the free energy to investigat e wh ich kinds of glassy phases can occur in 



site-disordered systems. The expansion is derived using results of section II 1. Since it includes only the lowest order 
overlap (cr Q cr/3), it has the same form as the Ginzburg-Landau expansion for the SK-model The difference lies in 
the prefactors, which now depend on the concentration of spins and on temperature. We show that these prefactors 
can change sign, giving rise to phase-transitions between different glassy phase. The phase diagram of the model is 
studied as a function of the prefactors. 



15 



In order to obtain the Ginzburg-Landau potential, we expand the free energy eq. (^) in the off-diagonal elements 
of p a p and q a /3, which vanish at the transition. In the following, we redefine q aa = p aa = and keep track of the 
diagonal elements pd and qd explicitly. We obtain 



pK = I d<ifcin(i " c ^ k ^ - lit j' j : pj - x> j )« 



~^2qal3Pa 



a/3 



n 

+ -^qdPd 
/ 



i / l 



+n log 2 + H 1 ~ tr s } E J\ ( 2 E P*P a <* a P ) ex P 



\ 



3=0 ■ 



a/3 



7^E CT ' 



where Jj is the j-th moment of the effective coupling (E 



d d k 



(27T) 



d^y 



(95) 



(96) 



As is discussed in section II 1 , this coupling includes the effects of clustering of spins which become important at 
higher spin densities. The saddle point equations for p a p and q a p are readily obtained from eq.(|95|). 



Pa0 = j2^ c ^ lj j(i^ 1 )^ 



3=2 



q al 3 = ^7/tri°cr Q cr^expX W = {<r a (Tp). (97) 
I 

(98) 

Using these equations, the p Q( a are eliminated from the free energy. We subtract the lowest order contribution, which 
gives the paramagnetic background, and rescale cf3 2 J2q a /3 — * qap- eq.(95j) now reads 



3=2 



1\ T J Jj 



where 



/=1 



v 



j=0 



J'! 2 



a/3 



We expand the free energy to fourth order. gives the following second-order contribution 
^E 7/tr s° E fci/Sife^crd^^^cr^exp \pdY a o 



(99) 



(100) 



(101) 



i ai/3ice 2 /32 L a 

The trace only has a non- zero value if the replica-indices are pair-wise equal. Because q aa = 0, ot\ = (3\ and ai = /?2 
give no contribution. Eq.(lOl) thus gives 



a 

where we have introduced the Z-spin moments Hkx—k m 

Wfei-k™ = ET' tr * )<7 i 1 ••• (7 m m exp 
I 

ECO (0 ( 
H TO fel • • • m fc m ( TO o J 

i 

(0 (0 



(102) 



(0 CO ' 



(103) 
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with mj? = tr a a k exp[p d <7 2 /2]. 

The third-order contribution from j — 3 in eq.(lOO) is derived in the same way. The replica- indices can be paired 
up in eight ways, each yielding a term X^? 3 )""- F rom 3 = 2 comes a term of the same form. 

Both j = 2 and j = 3 give the fourth order term ^2 aa {q 4 )aa, 3 ' = 4 gives: 

E 9a/3 8 ways. 

a/3 

^qlffVlj 48 ways, 

a/37 

E 1a00' r q 1 sqsa\ a¥:Til3 jL 5 48 ways, 

E tiM^ww - 12 ws (104) 

a/37<5 



Since E^^L) is of order n , the last term in eq. (104) will not contribute in the limit n — > 0. We implement 



the constraints on the indices by Kronecker delta's and then sum them out. For the third term in eq. (104), this for 
instance gives 

E 9apQPfQt5QSa\ a ^: Tt i3^:S = E QaplpyQjSqSaO- ~ &aey)0- ~ <W 

= E Vapqp-tq-fsqsa - 2 E vlpili + E q U- ( 105 ) 

af3~{6 a/3~{ a/3 

Combining these results with the contributions from the first term in eq. (p9|), we finally obtain the following free 
energy, which has the same form as for the SK model. 

w« - -~ E(« 2 )- - f E^ 3 u - 1 E & - 1 E - 1 E(? 4 )- ( 106 ) 

a a a/3 a/37 Q 

The prefactors in this expansion are functions of concentration and temperature. The prefactor of the quadratic 
term, r = (^22 ~ T 2 / CJ2) /2, vanishes at the spin glass temperature T g (c) = \fcJ2l^n.- The other prefactors in the 
free energy are given by 

w = ^222 - 2 T + 3/i22^y^, (107) 
3 1 

2/1 = ^^2222 + - M422, (108) 

V2 = fH22 - 3^2222, (109) 

rV 4 r 2 (j| + 2J2J 4 ) tj 3 

yj- + <^22 ^4 1- 4/1222—^- 

'2 J 2 J 2 



2/3=^2222 - 3 ^+2^22 ^ h 4/i 2 22 ■ (HO) 



The paramagnetic behavior is coded in the parameters pd and that satisfy the coupled mean field equations 
Pd = PJi and qd = P2- All information on clustering is contained in r, w, and the y's, so in the /i's and the Jj. In the 
limit c — > the /it's go to unity and for T ~ y^c the J3 and J4 terms vanish, so that one recovers the Ginzburg-Landau 
free energy of the SK-model. The important factors then are w = 1, y\ = 2/3, while the values of 2/2 (= —2) and 2/3 
(= 1) are irrelevant. When following the transition line T — T g (c) in the c — T phase diagram as function of c, it is 
seen that the higher /i's are rapidly oscillating functions. For instance, if J3 and J4 are much smaller than J2, then 
yi changes sign at c = 2.7% and at c = 4.3%, while w becomes negative at 6.7%. 

Based on these observations we are led to assume that the relevant physics near the phase transition(s) is still 



contained in the GL free energy eq. (10£). However, there is no reason to assume that w and y\ will always be 
positive. (A sign change of y\ occurs also in a Potts glass. j35|) Given the type of the lattice and the values of the 
spin-spin couplings, the c — T phase diagram may exhibit a limited number of special points (c*,T*) where either w or 
yi vanishes, and new phase boundaries originate. In fig. [I], a fictitious phase diagram is depicted for a system where 
w changes sign at some (c*, T*). 

In the limit n — ► 0, the free energy can be expressed in terms of the Parisi order parameter function q(x) in the 
usual way M. 
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PF = Si dx {l q2{x) - j«w t w + f - ^ir^) jf + f t2 ^} > ( m ) 

where 

T{x) = ^xq 2 (x)+q(x) dyq(y) + ^J^ dyq 2 (y). (112) 

The phasediagram of this model is determined by the parameters w and y\. The sign of the cubic prefactor 
determines whether the model exhibits a second or a first order transition. For positive y\ , the low temperature phase 
will be replica symmetry broken, for negative y\ replica symmetric. In the region where y\ is of order r 2 , a sixth 



order replica symmetry breaking term becomes relevant. We will discuss this case in section VII B. The parameters 
yi and j/3 are only important in the region where w <C 1, which we will discuss first. 

A. Behavior near first order type phase transitions 

Suppose we take w negative and y\ positive. As we lower temperature, the system undergoes a first order transition 
to a state which breaks replica symmetry. At positive w, on the other hand, we get a second order transition and the 
low temperature state will undoubtedly be different from the one at negative w (although it will also break replica 
symmetry, since we keep y% positive). We want to study the transition between these two low-temperature state s. We 



therefore take y± > and allow w to change sign. For w ~ ^fr <§; 1, q is also proportional to y/r, so eq. (Ill) then 



At positive w, we find the usual transition from the paramagnet (PM) to a Parisi-type spin glass phase (SG) with 
infinite replica symmetry breaking. At w < 0, however, the order parameter in the low-temperature phase only has 
one-step replica symmetry breaking (1RSB). Solutions of this form have previously been found by Gross, Kanter and 
Sompolinsky jj5| for the Potts glass and by Crisanti, Horner and Sommers (CHS) J3(| for a spherical spin glass model 
with p-spin interaction. We find a continuous transition between the SG and the 1RSB phase (fig. ^). 

The prefactors yi and 2/3, which are only important if w is small, must obey the following conditions: 2/3 < —y\ 
and y2 > max[0, —2yi — y^]. The first is related to keeping the plateau q\ positive, the second to doing this for the 
breakpoint x\ (q\ and x\ are defined below). 

1. The spin glass (SG I) solution 

First, we need to derive Parisi's infinite-RSB solution in the presence of the y2 and 2/3-term. We use the following 
ansatz for the order parameter function q(x): 

«<•' - { StfUo : 11131 

Because q{x) is continuous, it is completely determined by the saddle-point equation SF/Sq(x) = 0, 

rq(x) - wT(x) + ^-q 3 (x) - ^^q(x) f dyq 2 (y) 



+V3 

where 



T{x)R(x) + / dyT(y)q(y) + q(x) / dyT(y) 



= 0, (114) 



R{x) = xq[x)+ I dyq(y). (115) 



Both T(x) and R(x) are constant on the plateau at x > x\. 

The stability of the SG solution is determined by the sign of the "ergodon" -function A(x), which is an eigenvalue 
of the fluctuation-matrix 5 2 F[q]/Sq(x)5q(y). 



A(x) = r - wR(x) + ^q 2 (x) - ^±J^ J dyq 2 (y) + y 3 



R 2 {x)+xT{x) + / dyT{y) 



(116) 
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Differentiating the saddle-point equation eq. (114) with respect to x, we obtain q'(x)A(x) = 0. This implies that the 
SG solution is marginally stable on the interval < x < X\. Moreover, it is also marginally stable for x > x±, because 
q(x) is continuous at the breakpoint. This is of course a well-known feature of this state. 

The fourth derivative of the saddle-point equation with respect to x yields a second order differential equation for 
q nc (x). Its second derivative can be used to fix the height of the plateau q\. We obtain, 



Qnc{x) = — — ■ (117) 

3(?/i + 2/3^1) vfi + V3 X 



Finally, the breakpoint x\ is given by the implicit equation A(xi) = 0. 

) (xi - aatanh (^)) 

(118) 



Xl (3yi + 3y 3 - 2y 2 )xj + f 3 (yi+ yzx\) (x x - aatanh (^)) 



3(3/1 + 3/3x1) 18(yi + yzx{f 

where a = t/w 2 and a = \j —y\jy 3 . Near the transition, this yields x\ ~ 3?/icr. 

For a <C 1, w > 0, eq.(117) reproduces Parisi's solution of the SK model with x\ <C 1 (fig. ||a). When x\ approaches 



a, q(x) squeezes into a stepfunction, with the discontinuity located at x = x\ (fig.|2jf). This is the transition from 
infinite to one-step replica symmetry breaking. It occurs in the region w < 0, at a = <j c = l/6z/3 + yr>/18y 2 (l — a) 
(ng. |. 



2. One-step replica symmetry breaking solutions 

The order parameter in the 1RSB state has the form of a step function 

/ / ; < x < xi , . 

q{x) = 1 gi ; x x < x < 1 ( 119 ) 

Two conditions are needed in order to fix both the plateau q\ and the breakpoint X\. Next to the saddle-point equation 
dF/dqi = 0, we may impose either stationarity with respect to the breakpoint, dF/dxi — 0, or marginal stability of 
the ergodon on one of the plateaus, A /i = 0. We will consider both procedures. The first gives the proper static 
solution (or the solution for an infinitely slow cooling rate), while the second solution ought to occur in dynamics. 
The saddle-point equation dF/dqi = reads 

t - w { 1 - \x x J qi + \{yi- 2/ 2 (l - x x ) + 3y 3 (l - x±) + y 3 xj)ql = 0. (120) 



2 V 2 

The static solution must satisfy the additional equation 



2 \ 1 1 /3 3 \ 

1 - 3»i J ?i - iVill + 2^2(1 - xi)ql -y 3 {--x 1 + -x 2 1 jq 2 = 0. (121) 

The breakpoint of the marginally stable solution is fixed by 

3 1 

Ai = -r + wqi - -{y x + y 3 )q 2 + -y 2 (l ~ x x )q\ = 0. (122) 

In principle, the condition Ao = could also yield a dynamical solution, but this 1RSB state turns out to be unstable. 
We find the following static solution 

gf at - 3 +0 W \ 1 V (123) 
f 2/i + 32/ 3 xi (1 - ijClj 

This solution sets in from X\ = 1, at a negative value of a, cr^ tat = 1/9(3/1 + 3/3). The transition is first order in the 
sense that the order-parameter changes discontinuously and that the transition occurs before the AT instability of the 
PM phase. There is no latent heat, however. In fact, there is a problem, because the continuation of the PM phase 
has a lower free energy than the new 1RSB solution. The system nonetheless needs to go to the 1RSB state, because 
the PM becomes unstable at a — 0. This problem was also encountered by CHS. 
The ergodon mass of the static 1RSB is the same on both plateaus 
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A stat = Af at = + 2/3X 2 )(g stat ) 2_ (124) 

It goes unstable at the transition to the SG phase, cr = a c , where x\ becomes equal to a. 
The dynamic solution is given by 

„ mar s - TEl fi2K1 

91 - 22/1 + y 3 x 1 (3-x 1 )- (125j 

It also sets in a cci = 1, but at a higher temperature than the static 1RSB: cr" la ' s = 1/8(2/1 + 7/3). Close to cr" lalg , its 
free energy is lower than that of the PM (AF even has a finite slope, but there is no latent heat |37[] , |p8[ ). However, 
it becomes higher than -Fpm before the PM goes unstable at a = 0, so we have the same problem as with the static 
solution. 

The dynamic solution has Ai — 0, while the other eigenvalue is given by 

K^ = -\i.vt + y^l)(<ir s f- (126) 

It vanishes at the transition to the SG phase. 

The subject of dynamical transitions has made large progress since the start of this work, for a review see [Q. It 
is now understood that the marginality criterion signals dynamics in the highest, marginal states. This is the regime 
reached when first taking the thermodynamic limit and then letting time become large. In realistic systems one expects 
to go to lower states in the course of time, in a way that can be fixed, to some extent, by the cooling procedure. That 
happens in mean field spin glasses provided one considers dynamics at time scales exponential in the system size [ fl0"| . 
In other words, the marginality only points at the onset of a more interesting behavior. Dynamics in marginal states 
itself has a very limited meaning; it is a mean field artefact. The final result of these investigations has been the 
development of a picture for the thermodynamics of the non-equilibrium glassy state. Apart from the real tem per ature 
there occurs an effective temperature T e — T/xi, at which the system's slow modes are at quasi-equilibrium p(| 
This approach explains the old paradoxes related to the Ehrenfest relations and the Prigogine-Defay ratio J38| , that 
were the basis for the general belief that thermodynamics does not work for the glassy state. 



B. The transition from the Edwards- Anderson to the spin glass phase 

Another region of interest is where the prefactor of the replica symmetry breaking term, j/i, changes sign. For 
yi ~ —1, w ~ 1, the model is in the Edwards Anderson (EA) phase, which has a replica symmetric order parameter 
(fig.||e). At positive y±, it is in the SG phase. At the transition between these two phases, and higher order 

replica symmetry breaking terms become relevant. We therefore consider the following free energy 



^ = -\ E(? 2 )- - 1 E^ w - f E <£# - f E - 1 E ( 127 ) 

a/3 a/3 a/3 



The 2/2 and 2/3 terms are omitted, since they are not important for w ~ 1. We have included the most dangerous fifth 
and sixth order term. The full list of fifth order terms is 

J3 2T 5 J? 



Ef 5\ ( 1 , 1 rr 3^J2 + J3J4 . 1 r r 2^J2 + J 3 1 rj^Ji 2 

[Q Iota I YqA*22222 + ^1^22-L J5 r ^222± Jl h ^M2222-t -J2 ~ ^ 

E3 / 2\ 1 1 lrri^3 

Qafsil )a/3 I TM22222 - TJM4222 + T^/^442 + 



a0 



~ E <?a/30? 3 )aa ( -|/J22222 + ^4222 ) - E ^ ^ 3 ) 7 7 ^22222 (128) 



. 24' 

a/3 a/37 



The last of those is of order n 2 and does not contribute for n — > 0. The first and the third term are also neglected, 
since they hav e a s tructure similar to the cubic term, but are of higher order in r. It turns out that the j/5-term (the 



second of eq. (128)) can be absorbed into a redefinition of y\ by using the saddle point equation (q ) a p — ~2Tq a p/w. 
We therefore replace y\ by y\ = yi — 2ry^/w and omit the fifth order term. The prefactor of the most dangerous 
sixth order term is given by 
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15 15 1 15 1 1 

U& — — ^"M222222 + ^-^42222 — ^6222 — ^A*4422 + g^642 — 24Q^ 66 ' (129) 

We may expect new behavior when j/i and y§ become of the same order of magnitude. Writing j/i — 2y 4 r 2 /w 2 , we 
study the phasediagram as a function of 2/4 for r <C 1 and fixed. We find the following results: 

For ?/6 > 0, the model has a SG phase with nonzero lower plateau go- It interpolates between the EA phase, where 
go = 91, and the SG phase, where go = (figjj). 

For negative y§, there are two saddle points at the transition between the EA and the SG phase. The metastable 
saddle point consists of a 1RSB solution which sets in from the EA phase, followed by a transition to a new phase 
which we call SG/V (fig.||d). At larger j/4 there is a second transition from SG-ZVto SG. 

The order parameter of the stable saddle point is shown in fig.||c. This is also a new solution, which we call SGI/7. It 
has a lower free energy than the 1RSB and the SG/F solution. The SGIIIsets in from the EA phase with x\ = go — 0, 
and evolves into the SG at higher t/4, where go becomes equal to q±. 

The phasediagram for this region of parameter space is shown in fig. 

The solutions in this region are obtained in much the same way as t 
that there is no direct transition from the EA to the SG phase and then simply state the results for the phases which 
occur in between. 



re SG and the 1RSB. We therefore only show 



As n — > 0, the free energy eq. ( 127 ) takes the following form in terms of the order parameter function g(a;) 



0F = J X dx {^(x) - ~q{x)T{x) + - f 1 6 ^)) ■ (130) 

Taking a replica symmetric ansatz, we now find for the order parameter in the EA phase 

g E A = - + ^ £ r 4 + 0(r 5 ). (131) 

w w° 

The SG solution is again determined by the saddle point equation 5F/5q(x) = 0. From the second derivative of 
this equation, the non-constant part of q{x) is obtained 



1 

w \ ' "111 1 



x(q) = ^ ( 6y±q— - 20y 6 q 3 ) . (132) 



The upper plateau and the breakpoint of the SG order parameter are given by 

* = - + ^V + o(A 

w w° 

6y 4 - 20y 6 3 4 
xi=x{qi) = t° + 0(t). (133) 

From these equations, it is seen that there is no smooth transition from the EA to the SG phase for nonzero y§. The 
details of the various solutions which occur in the region \y$\ ~ \y $\ are listed below. For all solution of the SG-type, 



the non-constant part of the order-parameter is given by eq.(132) 
SG with go 7^ 

T , /*V 2^ 

w V 2y 6 



1x=qiSG 1 x 1 =x 1SG . (134) 



SG777 



1RSB 



<--i + iW i+ ^ +cM <i35) 

qi = -r+ %6g r 3 - 3 W g 2 r 2 ) - ^|(20g 3 r - 15wtft) + 0(r 5 ) (136) 



w 



3j/4 -4y 6 3 

x x = — r J . (137) 
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SGIV 



2w 2w V 2/6 

T . 2/4 - 2/6 , 4 3^ J/4 + 2y 6 2 2 2/6 /„ a Q 4\ , /n/ 5\ /i QO \ 

91 = - + 2 p— r +wq T a ) 3 q Q T (2<$t - 3wq Q + 0(t°) 138) 

xi = 2 v ± ^y 1 (r 3 + 2w?or2) _ ye (6 2r + g 3) + 0(r 4j 

C. Unsuccessful attempts 

Since the procedure for finding saddle points of the free energy is to guess an Ansatz and then to check if it works, 
we cannot guarantee that we have found everything there is. In particular, we have not found satisfactory solutions 
for the region w < and j/i < 0. Also, there may be extra solutions in the regions where we have already found 
stable saddle points, as the results for y\ <C 1 show. There are however some states which we definitely know not 
to be there, since we have unsuccessfully tried them. They are shown in fig. [| Most notably, we have not found a 
stable 2RSB state (fig. |a) for any values of the parameters. One-step and infinite replica symmetry breaking seem 
to be the only types of RSB that occur. Also 1RSB states with a non-zero lower plateau (in zero magnetic field) have 
been tried without success. States 6, c and d were thought to evolve from the 1RSB at small w, but none of them do. 
State b (SGI/) does occur in the Potts glass J35j. 




c, 

c 

FIG. 1. c — T phase diagram for a fictitious system with a line w(c, T) = 0. PM=paramagnet; SG=spin glass. 
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q(x) 
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SGIII 
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EA 
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q{x) 



SG/7 



d) 
q(x) 

qi 
qo 



f) 

q(x) 

qi 



A SGIV 



+ 



xo X\ 



1RSB 



X Xl X 

FIG. 2. The order-parameter function q(x) for various phases. The SG/7 solution has been found for the Potts glass by Gross 
em et. al., but does not occur here. 



r = r, 




FIG. 3. t ~ w phase diagram for a system with j/i > 0, yz < —y\ and j/2 > max[0, 
right to left it may appear in Fig. 1 around the point (c», T»). 



-2yi — y-i] ; with w increasing from 
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PM 



o yi 



EA J 


I SG 

g = 










2/4 = -2y 6 2/4 


= 



FIG. 4. j/i — r phase diagram for w > 0, 2/6 > 0. The function q(x) in the SG phase is drawn in figure 2. a for the case qo = 0. 
In the EA-phase q(x) is constant (no RSB). 



PM 



o w — » 




FIG. 5. 2/i — t phase diagram for w > 0, 2/6 < 0. In the SGIII phase q(x) is as in figure 2c. Dynamically this phase splits up 
in a 1RSB phase and a SG IV phase , see fig. 2b, d. 
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FIG. 6. Some unsuccessful attempts 
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VIII. SUMMARY 



We have presented the details of a new description for phase transitions in site-disordered magnets, a field that 
was started by us in 1992. As described in a letter of 1993 |To| , the approach starts from considering an infinity 
of macroscopic order parameters. This approach is different from, but probably equivalent to, the space-dependent 



order parameter field showing up usually, such as in the model of section IV . 

For the case of a random-site Ising spin glass, the central result is given by eqs. (|5l|), (|37|), (|l8|), and (|is|). 
Simplifications are found in the combined limit of small concentrations and large coordination number z. Here the 
random site problem generalizes the mean field model with random bonds, the spin glass phase of the Hopfield model 
and of the SK model. Only the pair overlap order parameter q a @ plays a role. In general, an infinity of multistate 

overlaps qal---a k , and conjugated overlaps Pa}---a k occurs. A resummation has been done, which is useful for low 
temperatures. This approach simplifies a bit for small concentrations, where a close relation is found to the spin glass 
model with diluted random bonds of Viana and Bray. In this regime the random site model exhibits a percolation 
threshold for short range systems and the "concentration scaling" T g ~ c for RKKY systems. 

We have generalized the approach to the case of vector spins. This has led to a criterion for the transition 
temperature at low concentrations. 

In dilute YGd single crystals, closely related to dilute AuFe and CuMn crystals, a thermodynamic transition to a 
state with "complete" spin density wave formation occurs |fi2|. The formation of incomplete spin density waves in 
metallic spin glasses have been observed by neutron scattering |12| . Therefore incomplete spin density waves have been 
presumed to be the cause of spinglass behavior in metals |43J . In our theory we indeed see that the effective coupling 
Jcs(k) = J(k)/(1 — c[3J{k)(q<i — qEA)) enhances the quantitative contribution of maxima in J(k). Nevertheless, 
thermodynamics involves all wavevectors (see, e.g. ([23])) without decisive role for the special wavevectors connected 
to spin density waves. As expected already by many, and confirmed by the experiments of Weissman incomplete 
spin density waves do not play a distinctive role in the spin glass phase. Although our present theory should be capable 
to explain why dilute YGd has spin density waves, while AuFe and CuMn have a spin glass phase, it is not clear to 
us how this should be shown. Let us recall that in the second order cumulant expansion (Gaussian approximation) 
the spin glass transition always precedes spin density waves or ferromagnetism ]l0| . 

Our present theoretical analysis also applies to the ill understood cluster glass or mictomagnetic phase. In metallic 
systems this is the phase observed between the low concentration spinglass and high concentration (anti-) ferromagnet, 
see Jl[] for a review. We have considered the Ginzburg-Landau expansion for our site-disordered field theory, within 
the second order cumulant expansion (Gaussian approximation in replica space). These expressions show that the 
prefactors of the cubic and quartic terms in the Ginzburg-Landau expansion may go through zero at certain points. 
The physically interesting cases are when either the prefactor of the cubic term vanishes, or the one of the quartic 
term that is responsible for replica symmetry breaking. We show that in the first case there is a transition to a phase 
with one step of replica symmetry breaking, while in the latter there may occur new spin glass phases between the 
standard spin glass and the replica symmetric Edwards- Anderson phase. 

It would be interesting to analyze the clusterglass experimentally in more detail, and in particular to look for new 
phases. Theoretically it is worth to analyze loop effects for the new phases. They may shed insight on the question 
whether fluctuations (de-) stabilize the new phases, though their effect may also be a mere shift of the transition 



point. This problem is partly related to the notorious problem of the renormalization of the AT-line. 44 Finally the 
question of reentrant phases should also fall within the scope of the theory presented here. 
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